Multiorder coherent Raman scattering of a quantum probe field 
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We study the multiorder coherent Raman scattering of a quantum probe field in a far-off-resonance 
medium with a prepared coherence. Under the conditions of negligible dispersion and limited band- 
width, we derive a Bessel-function solution for the sideband field operators. We analytically and 
numerically calculate various quantum statistical characteristics of the sideband fields. We show 
that the multiorder coherent Raman process can replicate the statistical properties of a single-mode 
quantum probe field into a broad comb of generated Raman sidebands. We also study the mixing 
and modulation of photon statistical properties in the case of two-mode input. We show that the 
prepared Raman coherence and the medium length can be used as control parameters to switch a 
sideband field from one type of photon statistics to another type, or from a non-squeezed state to a 
squeezed state and vice versa. 

PACS numbers: 42.50.Gy, 42.50.Dv, 42.65. Dr, 42.65.Ky 



INTRODUCTION 



The parametric beating of a weak probe field with a 
prepared Raman coherence in a far-off-resonance medium 
has been extensively studied Q, 0, 0, H It has been 
demonstrated that multimode laser radiation and in- 
coherent fluorescent light Q can be replicated into Ra- 
man sidebands. Since a substantial molecular coher- 
ence can be produced by the two-color adiabatic Raman 
pumping method the quantum conversion effi- 

ciency of the parametric beating technique can be main- 
tained high even for weak light with less than one photon 
per wave packet To describe the statistical proper- 
ties of a weak quantum probe and its first-order Stokes 
and anti-Stokes sidebands in the parametric beating pro- 
cess, a simplified quantum treatment has recently been 
performed 0. It has been shown that the statistical 
properties of the quantum probe can be replicated into 
the two sidebands nearest to the input line, in agreement 
with the experimental observations 0, • 

However, many experiments have reported the obser- 
vations of ultrabroad Raman spectra with a large num- 
ber of sidebands 0,0, [3 El In the experiments with 
solid hydrogen 0,13, at least two anti-Stokes sidebands 
and two Stokes sidebands have been observed. In the 
experiment with molecular deuterium 0, a large Ra- 
man coherence \p a b\ — 0.33 and about 20 Raman side- 
bands, covering a wide spectral range from near infrared 
through vacuum ultraviolet, have been generated. In 
rare-earth doped dielectrics with low Raman frequency 
and long-lived spin coherence, a substantial Raman co- 
herence \p a b\ — 0.25 and an extremely large number of 
sidebands (about 10 4 ) can also be generated [T(j- Broad 
combs of Raman sidebands 0- [J [J 0- El have been in- 
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tensively studied because they may synt hesize to subfem- 
tosecond @,H>E10] and subcycle [u| pulses. The gen- 
eration of broad combs of Raman sidebands has always 
been examined as a semiclassical problem. While classi- 
cal treatments are sufficient for many purposes, a quan- 
tum treatment is required when the statistical properties 
of the radiation fields are important. On the other hand, 
broad combs of Raman sidebands with similar nonclas- 
sical properties and different frequencies may find use- 
ful applications for high-performance optical communica- 
tion. Therefore, it is intriguing to examine the quantum 
aspects of high-order coherent Raman processes. 

In this paper, we extend the treatment of Ref. Q 
to study various quantum properties of multiorder side- 
bands generated by the beating of a quantum probe field 
with a prepared Raman coherence in a far-off-resonance 
medium. Under the conditions of negligible dispersion 
and limited bandwidth, we derive a Bessel-function solu- 
tion for the sideband field operators. We analytically and 
numerically calculate various quantum statistical charac- 
teristics of the sideband fields generated from a single- 
mode quantum input. We show that, with increasing 
the effective medium length or the Raman sideband or- 
der, the autocorrelation functions, cross-correlation func- 
tions, photon- number distributions, and squeezing fac- 
tors undergo oscillations governed by the Bessel func- 
tions. Meanwhile, the normalized autocorrelation func- 
tions and normalized squeezing factors of the single-mode 
probe field are not altered and can be replicated into a 
broad comb of generated multiorder Raman sidebands. 
We study the mixing and modulation of photon statisti- 
cal properties in the case of two-mode input. We show 
that the prepared Raman coherence and the medium 
length can be used as control parameters to switch a side- 
band field from one type of photon statistics to another 
type, or from a non-squeezed state to a squeezed state 
and vice versa. We also discuss two-photon interference 
in coherent Raman scattering. Although the multiorder 
coherent Raman scattering can produce a broad comb of 
sideband fields with different frequencies, it behaves in 
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many aspects as a beam splitter [1J, [15|, LL21 UiJ LU] • 

Therefore, in this paper, we also make comparison of 
this conventional device with our system as and when it 
is possible. 

Before we proceed, we note that, in related problems, 
the generation of correlated photons using the y^ 2 ^ and 

,(3) 
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y} A > parametric processes has been studied 
The correlations between the Stokes and anti-Stokes side- 
bands and the possibility of transferring a quantum state 
of light from one carrier frequency to another carrier fre- 
quency (multiplexing) have been discussed for resonant 
systems |2lj . 

The paper is organized as follows. In Sec. [H] we de- 
scribe the model and present the basic equations. In Sec. 
IIII we study various quantum characteristics of the side- 
band fields generated from a single- mode quantum input. 
In Sec. lIVI we discuss the quantum properties of the side- 
band fields generated from a two-mode quantum input. 
Finally, we present the conclusions in Sec. IVl 
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FIG. 1: Principle of the technique: Two classical laser fields 
drive a Raman transition of molecules in a far-off-resonance 
medium. The beating of a weak quantum probe field with the 
prepared Raman coherence produces new sideband fields. 



II. MODEL 

We consider a far-off-resonance Raman medium shown 
schematically in Fig.^ Level a with energy Lu a is coupled 
to level b with energy ujb by a Raman transition via inter- 
mediate levels that are not shown in the figure. We send a 
pair of long, strong, classical laser fields, with carrier fre- 
quencies u>_l and Wq , and a short, weak, quantum probe 

field Ei n , with one or several carrier frequencies, through 
the Raman medium, along the z direction. The tim- 
ing and alignment of these fields are such that they sub- 
stantially overlap with each other during the interaction 
process. The driving laser fields are tuned close to the 
Raman transition a <-> 6, with a small finite two-photon 
detuning 5, but are far detuned from the upper electronic 
states j of the molecules. We assume that all the fre- 
quency components of the input probe field are separated 
by integer multiples of the Raman modulation frequency 
uj m = LL>b — uj a — S. The drivingfields adiabatically pro- 
duce a Raman coherence p a b 0, 01 ■ When the probe field 
propagates through the medium, it beats with the pre- 
pared Raman coherence. Since the probe field is weak 
and short compared to the driving fields, the medium 
state and the driving fields do not change substantially 
during this step. The beating of the probe field with the 
prepared Raman coherence leads to the generation of new 
sidebands in the total output field -E ut . The frequencies 
of the sideband fields E q are given by u> q — u>q + qu) m , 
where q is integer and ujq is a carrier frequency of the in- 
put probe field. The range of q should be appropriate so 
that ijj q is positive. The probe field is taken to be not too 
short so that the Fourier-transformation limited broad- 
ening is negligible. We assume that the prepared Raman 
coherence p a b is substantial so that the spontaneous Ra- 
man process is negligible compared to the stimulated and 
parametric processes. Consequently, the quantum noise 
can be neglected. Unlike Ref. our model does not re- 



quire any restriction on the magnitude of the coherence 
as all Raman sidebands are included. When we take the 
propagation equation for the classical Raman sidebands 
|3, la] and replace the field amplitudes by the quantum 
operators, we obtain 

OE OE 

~dz + ~~c&t = i M u 9 E 9 + dq-iPbaE q -i + d q p ab E q+1 ). 

Here, u q and d q are the dispersion and coupling con- 
stants, respectively. We have denoted j3 q — Affiuj q /eoc, 
where Af is the molecular number density. 

We take all the sidebands to be sufficiently far from 
resonance that the dispersion of the medium is negligi- 
ble. In this case, we have u q = uo and d q = do. We 
write p ab = p exp[i(<f> - f3 m u z)}, where p a = \p ab \ and 
fim = N"hu> m /^o c i an d assume that po and </>o are con- 
stant in time and space. We change the variables by E q — 
£ q exp[i(/3 q uoz — q<f>o)]- Using photon operators, we can 
write £ q (z,t) = (2fiu q /e LA)V 2 J2 K b q (K,t)e iK ( z - ct l 
Here, L is the quantization length taken to be equal to 
the medium length, A is the quantization transverse area 
taken to be equal to the beam area, K is a Bloch wave 
vector, and b q (K, t) and bt(K, t) are the annihilation and 
creation operators for photons in the spectral mode q 
and the spatial mode K . For simplicity, we restrict our 
discussion to the case where each sideband field contains 
only a single spatial mode (with, e.g., K = 0). Then, 
Eq. yields 

db 

— = iiggbq^x + g q+1 b q+1 ), (2) 

where g q — (Afh j eo) y/Lo q u> q -i dgpo- For the medium 
length L, the evolution time is t = L/c. It follows from 
Eq. (J2J that the total photon number is conserved in time. 
Note that Eq. J2J) represents the Heisenberg equation for 
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the fields that are coupled to each other by the effective 
interaction Hamiltonian 

H=-n^g q+1 (b q bl +1 + b (!+1 bl). (3) 

9 

The interaction between the sideband fields via the pre- 
pared Raman coherence is analogous to the interaction 
between the transmitted and reflected fields from a con- 
ventional beam splitter |l9j . However, the two mecha- 
nisms are very different in physical nature. The most 
important difference between them is that the two fields 
from the conventional beam splitter have the same fre- 
quency while the sideband fields in the Raman scheme 
have different frequencies. In addition, the model Hamil- 
tonian @ involves an infinitely large number of Raman 
sidebands, separated by integer multiples of the Raman 
modulation frequency uj m . Despite these differences, the 
model © can be called the multiorder Raman beam 
splitter. The parameters g q t = g q L/c determine the 
transmission and scattering coefficients for the fields at 
the Raman beam splitter. 

We assume that the bandwidth of the generated Ra- 
man spectrum is small compared to the characteristic 
probe frequency loq. In this case, the g-dependence of 
the coupling parameters g q can be neglected, that is, we 
have g q = go = A/Tk^o^oPoAo- With this assumption, we 
find the following solution to Eq. @: 

b q (t) = Y,e l{q - q ' )7:/2 J q -A2g t)b q ,(0). (4) 
q' 

Here, Jk is the fcth-order Bessel function. The expression 
Q for the output field operators is a generalization of the 
Bessel-function solution obtained earlier for the classical 
fields 0, ■ The number of generated Raman sidebands 
is characterized by the effective interaction time 2g$t or, 
equivalently, the effective medium length kL, where 

2.9o 2fi 

k= = — Nuj doPo- (5) 

c e c 

The coefficient k characterizes the strength of the para- 
metric coupling and is proportional to the prepared Ra- 
man coherence po, that is, to the intensities of the driv- 
ing laser fields. The Bessel functions Jk{nL) are the 
transmission (fc = 0) and scattering (fc ^ 0) coeffi- 
cients for the Raman sidebands, similar to the transmis- 
sion and reflection coefficients of a conventional beam 
splitter. The assumption of limited bandwidth requires 
kL < uJo/uim, that is, (2h/eoc)Afu m d p L <C 1 0, @. 
In what follows we use the expression Q to calculate 
various quantum statistical characteristics of the side- 
band fields, namely, the autocorrelation functions, the 
two-mode cross-correlation functions, the squeezing fac- 
tors, and the relation between the P representations of 
the output and input states. 



III. SINGLE-MODE QUANTUM INPUT 

In this section, we consider the case where the input 
probe field has a single carrier frequency ujq. In other 
words, we assume that the sideband q = is initially 
prepared in a quantum state and the other sidebands 
are initially in the vacuum state. The density matrix of 
the initial state of the fields is given by 

Pin = 1^11 ( 6 ) 

A. Autocorrelation functions 

We study the autocorrelations of photons in the gen- 
erated Raman sidebands. We use Eq. Q and apply the 
initial density matrix JBJ to calculate the normally or- 
dered moments (b\ n bq) of the photon-number operators 
n q = Wqbq. The result is 

(6j B Sj) = Jf(«L)(SS"(0)S5(0)). (7) 

In particular, the mean photon numbers of the sidebands 
are given by 

(n q ) = J*{KL)(n in ). (8) 

Here, n;„ = Sg(0)&o(0) is the photon-number operator for 
the input field. The nth-order autocorrelation functions 
of the sidebands are defined by T q n) = (b\ n b r q L ) - {b\b q ) n . 
From Eqs. and ©, we find 

T^=J^{KL)Tt\ (9) 

where = (6J n (0)6ft(0)) - (6J(0)S (0)) n is the nth- 
order autocorrelation function of the input field. 

Equations © and @ indicate that, when we increase 
the effective medium length kL or the sideband order q, 
the mean photon number (n q ) and the autocorrelation 

(n) 

function T q undergo oscillations as described by even 
powers of the Bessel function J 9 (kL). Such oscillatory 
behavior is illustrated in Fig. [21 When the sideband order 
q is higher, the onset of (n q ) occurs later [see Fig. EJa)] 

and, hence, so does the onset of T q n ^ [see Fig. Etc)]. For 

a fixed q, both (h q ) and reach their largest values 
at the same optimal medium length L q — x q /n, where 
x q is the position of the first peak of J q (x). The higher 
the sideband order q, the larger is the optimal length L q 
and the smaller are the maximal output values of (h q ) 
and T q n) [see Figs.Eta) and Etc)]. Figures Etb) andEtd) 

show that {h q ) and T q n ^ are substantially different from 
zero only in the region where |gj is not too large compared 

to kL. For a given kL, both (n q ) and achieve their 
maximal values at q w kL. 

The normalized nth-order autocorrelation functions of 
the sidebands are defined by g^ = (b^ l b™) / (b^bq) 71 . 
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FIG. 2: 



Mean photon number (h q ) (first row) and second- 
order autocorrelation function P, 2 ' (second row), both scaled 
to their initial values for the probe field, as functions of the 
effective medium length kL (first column) and the sideband 
order q (second column). In (a) and (c), the sideband order is 
(solid line), 1 (dashed line), and 5 (dotted line). In (b) and 
(d), the effective medium length is kL = 5. In (c), we have 
amplified r£ 2) / r ta (dashed line) and rf } /F$ (dotted line) 
by 5 and 10 times, respectively. In (b) and (d), the negative 
side of the q axis is not shown because the functions plotted 
are symmetric in q. 



These functions characterize the overall statistical prop- 
erties, such as sub-Poissonian, Poissonian, or super- 
Poissonian photon statistics, regardless of the mean pho- 
ton number. Unlike the mean photon number (n q ) and 
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the normalized au- 



the autocorrelation function r 

(n) 

tocorrelation function g q does not oscillate when we 
change the effective medium length kL or the sideband 
order q. Indeed, with the help of Eq. (JJJ, we find 



where g ; 



(n) 



<^ n (0W))/(^(0)6o(0)r 



>) = n (n) 
Hq y m ' 



(10) 



Equation i|10[l indicates that the generated sideband 
fields and the probe field have the same normalized auto- 
correlation functions, which are independent of the evo- 
lution time and are solely determined by the statistical 
properties of the input field. In other words, the nor- 
malized autocorrelation functions of the probe field do 
not change during the parametric beating process and 
are precisely replicated into the comb of generated side- 
bands. Such a replication of the normalized autocorrela- 
tion characteristics can be called autocorrelation multi- 
plexing. In particular, if the photon statistics of the input 
field is sub-Poissonian, Poissonian, or super-Poissonian, 
the photon statistics of each of the sideband fields will 
also be sub-Poissonian, Poissonian, or super-Poissonian, 
respectively. This result is in agreement with the experi- 
ments on replication of multimode laser radiation and 



broadband incoherent light |3| ■ The ability of the Raman 
medium to multiplex the autocorrelation characteristics 
is similar to the ability of a conventional beam splitter 

G3. 

It is not surprising that the normalized autocorrela- 
tion functions of the probe field are replicated into the 
sidebands in the parametric beating process. Such a 
replication is possible because the medium is far off res- 
onance and the quantum probe field is weak compared 
to the driving fields. Under these two conditions, the 
photon annihilation operators are linearly transformed 
as described by the linear differential equation This 
equation shows that the photon annihilation operators 
are not mixed up with the creation operators, and there- 
fore the evolution of the annihilation operators is linear 
with respect to the initial annihilation operators. 

We should, however, emphasize here that the repli- 
cation of all of the normalized autocorrelation functions 
9in* (f° r ai l orders n) of the input field does not mean the 
replication of the quantum state /5j° . In fact, the oscilla- 
tions of the mean photon number (n q ) and the autocor- 

(n) 

relation function T q indicate that the photon-number 
distributions and consequently the quantum states of the 
sidebands evolve in a rather complicated way and are 
quite different from those for the input field. A sepa- 
rable state at the input can produce an entangled state 
0,|23. In addition, cross-correlations between the side- 
bands can be generated from initially uncorrelated fields, 
and a Fock state at the input does not produce sideband 
fields in isolated Fock states [see below]. 



B. Cross-correlation functions 

We study the correlations between the generated Ra- 
man sidebands. For two different sidebands k and I 
(k 7^ /), we have 



(h k fii) = J%(KL)J?(KL)(n in (n in - 1)). 



(11) 



The cross-correlation function for the two sidebands is 
(fikfii) - {h k ){hi). Using Eqs. $ and 



defined by rjy 



ifTTJl. we find 



.(2) _ 
kl ~ 



4(KL)Jf(KL)T\ 



(2) 



(12) 



When we extend T^? for k = I, we have = 1^ . 

_____ /<2~\ 

According to Eq. I|12|) . the cross-correlation function T kl 
oscillates when we change the effective medium length kL 
or the sideband orders k and I. Such oscillatory behavior 
is illustrated in Fig. |3J 

The normalized cross-correlation function is defined by 



9ki = (fi-kfii) I {{fik) (hi)) . Unlike the function T y kl ' , the 
(2) 

normalized function g kl does not oscillate. Indeed, we 
find the relation 



^(2) 



(2) 

9ki 



9 



(2) 
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Sideband order / 

(2) 

FIG. 3: Cross-correlation function T kl , scaled with the 
second-order autocorrelation function of the input field, as 
a function of (a) the effective medium length kL and (b) the 
sideband order I. Here, the sideband order k is fixed to 1. 
In (a), the sideband order I is (solid line), 2 (dashed line), 
and 5 (dotted line). In (b), the effective medium length is 
kL — 5. In (b), the negative side of the I axis is not shown 
because the function plotted is symmetric with respect to the 
sideband orders. 



As seen from the above relation, the normalized cross- 
correlation functions g kl for all possible sideband pairs 
{k, I) are equal to each other, to the normalized second- 

(2) 

order autocorrelation function g q for each sideband q, 
and to the normalized second-order autocorrelation func- 

(2) (2) 

tion <7 in of the input field. When g in ^ 1, that is, when 
the photon statistics of the input field is non-Poissonian, 

(2) 

we obtain g kl ' ^ 1, a signature of cross-correlations be- 
tween the sidebands. Such correlations are generated 
although the sidebands are initially not correlated. In 
particular, if the input field has a sub-Poissonian pho- 

(2) 

ton statistics (g in < 1), anti-correlations between the 

(2) 

sidebands (g kl < 1) will be generated. Note that the 
conventional beam splitters also have a similar property 
[l9|. The generation of cross-correlations between the 
sidebands indicates that the quantum states of the gen- 
erated sidebands are different from that of the input field. 
We emphasize that the anti-correlation generation can- 
not be explained by the classical statistics of the fields 



with positive P functions although the sideband dynam- 
ics is linear with respect to the field variables. 

C. Photon- number distributions 

To get deeper insight into the quantum properties of 
the generated Raman sidebands, we derive the photon- 
number distributions for the output fields. The joint 
photon-number distribution for the output fields is de- 
fined by Ps({n q }) = {{n q }\p ou t\{n q }) i where p out = 
U(L/c)pi n U'(L/c) is the density matrix of the output 
state. Here, U(L/c) = exp(—iHL/hc) is the evolution 
operator. Using Eqs. @ and JBJ, we find 

ft(K» =ft»Wrr^-rII j X kL )> ( 14 ) 

ll q n q- q 

where Pi n (n) = o{n\p^ \ n )o is the photon-number distri- 
bution of the input field, and N = ^ q n q is the total 
photon number. From Eq. (|14|) . the marginal photon- 
number distribution p q (n) for the sideband q is obtained 
as 

pM = E^r^ 1 - J >L)\ k PUn + k). 

Clearly, p q (n) is in general different from Pi n (n). Thus, 
despite the replication of the normalized autocorrelation 
functions, the photon-number distribution of the probe 
field is not replicated into the sidebands. 

We examine several particular cases. First, we consider 
the case where the probe field is initially prepared in 
a coherent state |a)o- This state is characterized by a 
Poisson distribution Pi n (n) = e~ N N n jn\ for the photon 
number, where jV = |a| 2 . In this case, Eqs. (|14fl and i|15|) 
yield P^({n q }) = Y[ q p q (n q ) and 

n™ 

p( n ) =e -ft._l > (16) 

respectively, where n q — NJ^(kL). Thus, the generated 
sidebands are not correlated, and the marginal photon- 
number distributions for the individual sidebands remain 
Poisson distributions during the evolution process. 

Second, we consider the case where the probe field 
is initially in a Fock state \N} . In this case, we have 
Pin(n) = <5n,AT. Therefore, Eq. lfT3|l yields 

p q (n) = Jf {kL)[1 - J>L)] N - n nK *l n)l (17) 

for n < N, and p q (n) = for n > N. Meanwhile, 
Eq. JTIJ yields P^({n q }) ^ Y[ q Pq{n q ), that is, the 
joint photon-number distribution is not a product of the 
marginal photon-number distributions for the individual 
sidebands. Thus, the sidebands are correlated. They are 
not generated in isolated Fock states. 
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Finally, we consider the case where the probe field 
is initially in a thermal state, which is characterized 
by a Boltzmann photon- number distribution p- ln (n) = 
N n /{N + l) n+1 . In this case, Eq. {T5J yields 



(n q + 1) 



n+l ' 



(18) 



where n n 



u q — NJ 2 (kL). Thus, the marginal photon- 
number distributions for the individual sidebands remain 
Boltzmann distributions during the evolution process. 
However, according to Eq. (|14f> . we have P^({n q }) ^ 
Y\ q Pq{n q ), a signature of correlations between the gener- 
ated sidebands. 



D. Squeezing 

We examine the squeezing of the field quadratures. A 
field quadrature of the gth mode is defined by X q = 
ble^+bge-^. We say that the gth mode is in a squeezed 
state if there exists such a phase <p that ((AX q ) 2 ) < 1 
or, equivalently, S q < 0, where S q = ((AX q ) 2 ) - 1. The 
squeezing degree is measured by the quantity — S q . Note 
that the relation between the squeezing factor S q and the 
conventional squeezing parameter r q is S q = er 2r i — 1. 
In terms of the photon operators, we have 



S q = 2[(6t6 9 )-(St)(S 9 )] + [«S^)-<S g ) 2 )e- 2 ^+c.c.]. (19) 



Using Eqs. and ©, we find 

(S 9 ) = e^' 2 J q {nL)(b {Q)) 



and 



(b 2 q )=e^J 2 ( K L)(b 2 (0)). 



(20) 



(21) 



When we insert Eqs. ©, HO}, and (|3TJ into Eq. JTJJJ, we 

obtain 



S q {ip + QTr/2) = J 2 ( K L)S in (ip). 



(22) 



Here, S[ n (ip) denotes the squeezing factor for the <p- 
quadrature of the input field. Equation shows that, 
if S in (ip) < 0, then S q (ip + qTr/2) < 0. Thus, if the input 
field is in a squeezed state, then the generated sidebands 
are also in squeezed states. In other words, the squeez- 
ing of the input field is transferred to the comb of gen- 
erated sidebands during the parametric beating process. 
The squeezing factor S q ((f + ?7r/2) of the sideband q is 
reduced from the input squeezing factor Si n (ip) by the 
factor J 2 (kL). Unlike the case of linear directional cou- 
plers and beam splitters , the squeezing degree of the 
probe field cannot be completely transferred to the Ra- 
man sidebands. This difference is due to the fact that the 
linear directional coupler and the beam splitter involve 
only two output modes while the multiorder coherent Ra- 
man process involves many more output modes. Note 



that the phase of the squeezed quadrature of the side- 
band q changes by qTr/2. This means that the squeezed 
quadrature of a generated even-order (odd-order) Raman 
sideband is parallel (orthogonal) to that of the input field. 

We introduce the normalized squeezing factor s q — 
S q /(n q ). We find the relation 



s q {ip + qTr/2) = s- m ((p), 



(23) 



where s- m ((p) = S- in ((p) / (n- m ) is the normalized squeez- 
ing factor for the input field. Equation <2.''>l) implies 
that, besides a shift of the quadrature phase angle, the 
normalized squeezing factor for the input field is repli- 
cated into the comb of generated sidebands. This result 
can be used to convert squeezing to a new frequency, 
i.e., to perform squeezing multiplexing. The relation 
S q (ip + qir/2)/ Si n ((p) = (n g )/(n; n ) indicates that the kL- 
and (/-dependences of the squeezing factor S q (ip + qTr/2) 
are similar to those of the mean photon number (n q ) [see 
Figs. Eta) and|2{b)]. Note that, when the input field is in 
a coherent state, the sideband fields have no squeezing. 
This property is similar to the case of four-wave mixing 
but is unlike the case of degenerate parametric down- 
conversion, where perfect squeezing can in principle be 
obtained. Since squeezed states are nonclassical states, 
the ability of the Raman medium to multiplex squeezing 
from a probe field to its sidebands is a quantum property 
that cannot be described by the classical statistics of the 
fields with positive P functions. 



E. Quantum states of the output fields 

We calculate the quantum state of the output fields for 
several cases. First, we consider the case where the input 
sideband is initially in a coherent state |a)o- The state 
of the fields at the input is 

l*in) = Ho II l°>« = e- |Q|2/2 e^ (0) |0). (24) 

The state of the fields at the output is given by l^out) = 
U(L/c)\V in ). Since U(L/c)\0) - |0), we have |* out ) = 

e -\a\ 2 /2 e abl(-L/c)\ Q y Using Eq _ we fi nc J 



i*out)=iK(L/ C )})=ni a «( L / c ))^ 



(25) 



Here, \a q (L/c)) q is a coherent state of the qth mode, with 
the amplitude 



a q {L/c) = aJ q (KL)e lq7T/2 . 



(26) 



Thus, a probe field in a coherent state can produce side- 
band fields that are also in coherent states. Such a 
process can be called coherent-state multiplexing. This 
property of the Raman medium is similar to the case of 
conventional beam splitters [l9l | . 
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Second, we consider the case where the input sideband 
is initially prepared in a Fock state \N) . The state of 
the fields at the input is written as 

|* i n) = |JV)oni°)9 = ^o JV (0)|0). (27) 

The output state of the fields is given by l^out) = 
(7V!)- 1 / 2 SS A '(-i/c)|0). With the help of Eq. ©, we find 

|*out) = £ C S}IW>- (28) 
{««} 

Here. 

y 1 lq Q q 

for E 9 n q = N, and c[% } = for J2 q ™q ¥= N. When 
N 7^ and kL ^ 0, the output state (|28|) is, in general, 
a multipartite inseparable (entangled) state. 

In a particular case where the input state of the probe 
field is a single-photon state, i.e., N = 1, Eqs. (|28Jl and 
(P| yield 

|*out) =J2 eiq * /2j «(*L)\l q ). (30) 

Here, |l g ) is the quantum state of a single photon in the 
sideband q with no photons in the other sidebands. In 
this case, the entanglement between two different side- 
bands k and I at the output can be measured by the 
bipartite concurrence Cki = 2\J} s (kL)Ji(kL)\ > see [22j. 

Finally, we consider the case where the input sideband 
is initially in an incoherent mixed state 

/4 0) =I> n (n)(|n>H)o. (31) 

n 

With the help of Eq. 1)28(1 . the density matrix of the out- 
put state of the fields is found to be 

Pout- £ Ptn(N)cl% } Cl$\{n g }){{n' q }\. (32) 

N,{n q },{n' q } 

To obtain the reduced density matrix p^ t for an arbi- 
trary sideband g, we trace the total density matrix l|32|) 
over all sidebands except for the sideband q. Then, we 
find 

n 

where the marginal photon- number distribution p q {n) is 
given by Eq. (|15fl . As seen, the reduced state of each side- 
band is also an incoherent superposition of Fock states. 
However, if the input state is different from the vacuum 



state, then, for kL ^ 0, we have p out ^ T\ q Pouti a s ig~ 
nature of correlations between the generated sidebands. 
Moreover, the total density matrix l|32l) of the output 
fields contains nonzero off-diagonal matrix elements in 
the Fock-state basis. 

In a particular case where the initial state of the 
probe field is a thermal state, i.e., p^ = J2 n [ Nn /i N + 
l)" +1 ](|n)(n|)o, the reduced state of each generated side- 
band is also a thermal state, namely, 

^ = E m (!"><"!)«• (34) 

Here, n q — NJ 2 (nL) is the mean photon number for the 
sideband q. The reduced thermal states of the generated 
sidebands are, however, not isolated from each other. 

IV. TWO-MODE QUANTUM INPUT 

A far-off-resonance medium with a substantial Raman 
coherence, prepared by two strong driving fields, can ef- 
ficiently mix and modulate the quantum statistical prop- 
erties of the sideband fields. To understand this mecha- 
nism, we study the case where the input probe field has 
two carrier frequencies, ojq and uj u = oJo + vuj m , separated 
by an integer multiple v of the Raman modulation fre- 
quency uj m . We assume that the Raman sidebands and 
v are initially in independent quantum states /3;° and 

p\^\ respectively, while the other sidebands are initially 
in the vacuum state. The density matrix of the initial 
state of the fields is given by 

An = 4° ) ®4n ) ® II (l°X°IV (35) 

Here, v ^ 0. 

A. Modulation of photon statistics 

We study the mixing and modulation of photon statis- 
tics of the sideband fields. When we use Eq. Q to calcu- 
late the mean photon numbers of the sidebands generated 
from the initial state (|35() . we find 

(b q b q ) = J 2 q ( K L)(bl(0)b (0)) + J 2 q _ v ( K L){bi (0)6,(0)) 
+ J q (K,L)J q ^ v (nL) 

x [e-W2(St( ))(S 1/ (0))+c.c.]. (36) 
Furthermore, we find 
(b\ 2 b 2 q ) = 4(*L)<6t2( )&g(0)) + J^(nL){bl 2 (0)bl(0)) 
+ 4 J 2 {nL)J 2 q _ u {kL) (SJ (0)So(0)) & (0)6,(0)) 
+ [2e-^/ 2 J 3 q (KL)J q ^(KL)(bl 2 (0)b (Q))(K(0)) 
+ 2e— I 2 J q (kL) (kL) (bl (0)) (St (0)bl (0)) 
+ e— J 2 ( K L)J 2 q _J K L)(bl 2 (0))(bl(0)) +c.c.]. (37) 
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Hence, the second-order autocorrelation function T q 
(bl 2 b 2 ) — (b„b q ) 2 of the sideband q is found to be 



(2) 



2.0 



r(2) 

q 



= J^kL)T {2 \o) + J^L)T( 2 \0) 
+ 2J 2 (kL)J 2 _ v (kL) Aq 
+ [e-^J*( K L)J 2 q _„( K L)A 1 



+ 2e- l ^ /2 J q (KL)j3( K L)A 3 + ex.], 



(38) 



where 

A = 

At -- 

A 2 = 

A 3 = 



(6S(o)8o(Q))<8t(o)^(o)>-|(5o(o))| a KS v (o))| 2 , 
(SS 2 (o))(S 2 (o))-<SJ(o)) 2 <S !/ (o)) 2 , 

<S.(0))[(S o 2 (0)S o (0)> - (6t( )S (0)>(fo (0)>], 
<SJ(0))[(6t (0)& 2 (0)> - (St(Q)S,(0))(6,(0))]. (39) 



The first two terms on the right-hand sides of Eqs. (13611 . 
l|37jl. and Ij38(l are the contributions of the individual in- 
put sidebands and v. The other terms result from the 
interference between the two interaction channels. 

Unlike the case of single- mode input, in the case of 
two- mode input, the normalized second-order autocorre- 



lation function g q ' = 1 + T y q ' / (n q ) 2 depends, in gen- 
eral, on kL and q. Such behavior is illustrated in Fig. ^ 
When kL is such that J q {nL) = or J q - v {nL) — 0, we 

have g q 2 ^ = c/u\o) or g q 2 ' — g^\o), respectively. Con- 
sequently, if the two input sidebands have different nor- 

(2) (2) 

malized autocorrelation functions, i.e., gfi 

then, with increasing kL or q, the normalized autocor- 

(2) 

relation function g q will oscillate between the values 

<7q 2 ^(0) and g^(0) [see Fig.0]. In particular, if the pho- 
ton statistics of one of the input fields, e.g., the sideband 
0, is sub-Poissonian ' (0) < 1] and that of the other in- 

(2) 

put field is super-Poissonian [g v (0) > 1], then each gen- 
erated sideband q will have complex statistical properties 

(2) 

and will oscillate between sub-Poissonian [g q ' < 1] and 
super-Poissonian [g q 2 ^ > 1] photon statistics [see Fig. 
Using the prepared Raman coherence po or the medium 
length L as a control parameter, we can switch a side- 
band field from super-Poissonian photon statistics to sub- 
Poissonian or vice versa. Similar modulation of photon 
statistics has been demonstrated in a linear directional 
coupler |l8[ . 



(2) 



B. Modulation of squeezing 

We study the mixing and modulation of the squeezing 
properties of the sideband fields. When we use Eq. J3J to 
calculate the amplitudes (b q ) and (b 2 ) of the sidebands 
generated from the initial state (|35|l . we find the expres- 
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FIG. 4: Normalized second-order autocorrelation function 
g q 2 ^ as a function of (a) the effective medium length kL and 
(b) the sideband order q in the case of two-mode input. The 
input sideband is initially prepared in a Fock state, with 
5 photons. The input sideband 1 is initially prepared in a 
thermal state, with 1 photon in average. In (a), the side- 
band order is —1 (solid line) and 2 (dashed line). In (b), the 
effective medium length is kL — 5. 



(S g ) = e^l 2 J q ( K L)(b Q (0)) + e^-^/ 2 J q ^( K L)(K(0)) 

(40) 

and 

(S 2 ) = e^J 2 (/ t £)(6 2 (0)) + e i ^>J 2 _ I/ ( K L)<S 2 (0)) 

+ 2e i <«-"/ 2 > *J q (nL) J q - V (kL) (Sq (0)) (K (0)) . 

(41) 

We insert Eqs. JHSJ) , 613), and (@T} into Eq. JE). Then, 
we obtain the squeezing factor 

S^ + qn/2) = J 2 (KL)St\v) 

+ J 2 _„(kL)SW(v + vit/2), (42) 

where Sg and S„ arc the initial squeezing factors of 
the sidebands and v, respectively. As seen, the squeez- 
ing factor S q of the sideband q is a superposition of the 

input squeezing factors Sq" 1 ^ and Sv, taken with the 
quadrature phase shifts —qn/2 and — (q — u)tt/2, respec- 
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FIG. 5: Normalized squeezing factor s q (q-K /2) as a function 
of (a) the effective medium length kL and (b) the sideband 
order q in the case of two-mode input. The input sideband 
is initially prepared in a squeezed vacuum state, with the 
squeezing parameter r = 1 and the phase 8 — 0. The input 
sideband 1 is initially prepared in a thermal state, with the 
mean photon number 1. In (a), the sideband order is — 1 
(solid line) and 2 (dashed line). In (b), the effective medium 
length is kL — 5. 

tively, and weighted by the factors J 2 (kL) and J 2 _ v (nL), 
respectively. 

Unlike the case of single- mode input, in the case of 
two- mode input, the normalized squeezing factor s q (ip + 
qn/2) = S q (ip + qn/2)/(h q ) varies, in general, with kL 
and q. Such behavior is illustrated in Fig. [S] When 
kL is such that J q (nL) = or J q - u (nL) = 0, we have 

s q (cp+qw/2) = s$,(<p+uir/2) or s q (ip + qir/2) = sfi n \cp), 
respectively. Consequently if the normalized squeez- 
ing factors Sq (y) and s„ n \(p + vit/2) of the two in- 
put fields are different, the normalized squeezing factor 
s q (ip + q7r/2) will oscillate between the values SQ n \tp) 
and Sv n \ip + b l ir/2). In particular, if one of the two input 
fields, e.g., the sideband 0, is squeezed [sq (ipo) < 0] and 
the other input field is not squeezed [s„ n \(p) > 0], then, 
each generated sideband q will have complex squeezing 
properties and will oscillate between a squeezed state 
[s q (ipo +qn/2) < 0] and a non-squeezed state [s q (<p) > 0], 
see Fig. \5\ Using the prepared Raman coherence po or the 
medium length L as a control parameter, we can switch 



FIG. 6: (a) Squeezing factor S q {qn/2) and (b) mean photon 
number (n q ) as functions of the effective medium length kL 
in the case where the sideband is initially prepared in a 
weak squeezed vacuum state and the sideband 1 is initially 
prepared in a strong coherent state. The parameters for the 
initial states of the input sidebands are r = 1, 9 = 0, and 
a = 20. The curves are calculated for the sidebands (solid 
line), 1 (dashed line), and 3 (dotted line). 



a sideband field from a non-squeezed state to a squeezed 
state or vice versa. Note that a similar result has been 
obtained for a linear directional coupler [l8|. 

We analyze a particular case where the sideband is 
initially in a squeezed vacuum state |£)o = cxp[(£*&Q — 
£&J 2 )/2]|0)o and the sideband v is initially in a coherent 
state \a) u . Here, £ = re %e is a complex number, the 
modulus r = |£| characterizes the amount of squeezing, 
and the phase angle 9 characterizes the alignment of the 
squeezed vacuum state in phase space. For the input 
sideband 0, we have [H (S£(0)S o (0)) = sinh 2 r, (6 (0)) = 
0, and SQ in \(p) = cosh 2r— cos (2ip—9) sinh 2r— 1. For the 
input sideband u, we have (6^(0)6 l/ (0)) = \a\ 2 , (b v (Q)) = 
a, and S^ in) (if) = 0. Then, we find from Eq. that 
the mean photon number of an arbitrary sideband q is 

{b\b q ) = J q 2 ( K L)smh 2 r + J 2 q _ v {nL)\a\ 2 . (43) 

We find from Eq. I|42(l that the maximal squeezing of 
the sideband q occurs in the </? g -quadrature where tp q = 
8/2 + qn/2. The corresponding value of the squeezing 
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factor is 

S q (<p g ) = J 2 q ( K L)(e- 2r -1). (44) 

As seen from Eq. (|44|) , squeezing can be transferred from 
the initial squeezed vacuum state of the sideband to 
the other sidebands. The squeezing factors of the side- 
bands are independent of the amplitude a of the initial 
coherent state of the sideband v. Meanwhile, the mean 
photon number of each sideband is governed not only by 
the squeezing parameter r of the initial state of the side- 
band but also by the amplitude a of the initial state of 
the sideband v. Using this fact, we can manipulate to get 
optimized mean photon numbers and squeezing degrees 
of the sideband fields at the output as per requirement. 
In particular, we can convert squeezing from a weak field 
to a much stronger field. To illustrate this possibility, we 
plot in Fig.|S|the squeezing factor S q (qn/2) and the mean 
photon number (n q ) as functions of the effective medium 
length kL for the parameters r = 1, = 0, v = 1, and 
a = 20. In this case, the most negative value of the input 
squeezing factor Sq (tp) is achieved at cp = and is given 

by SQ ln \o) = e~ 2r — 1 = —0.86, indicating the squeez- 
ing degree 86%. The mean photon number of the input 
squeezed vacuum state is (6q(0)6o(0)) = sinh 2 r = 1.38, 
rather small. The solid lines in Fig. [S] show that the 
sideband 0, initially prepared in a weak squeezed vac- 
uum state, can be significantly enhanced while keeping 
its squeezing degree substantial. Meanwhile, the dashed 
lines show that, for kL = 1.84, the sideband 1, initially 
prepared in a strong coherent state, is squeezed by about 
29% and has the mean photon number of about 41. Sim- 
ilarly, the dotted lines show that, for kL = 4.2, a gen- 
erated new sideband 3 is squeezed by about 16% and 
has the mean photon number of about 39. Thus, from a 
weak squeezed field at the input, we can obtain other out- 
put squeezed fields that have smaller but still substantial 
squeezing degrees, much larger mean photon numbers, 
and different frequencies. 



C. Two-photon interference 

We show the possibility of quantum interference be- 
tween the probability amplitudes for a pair of photons 
with different frequencies in the coherent Raman process. 
We assume that the sidebands and 1 are initially pre- 
pared in independent single-photon states. This initial 
condition corresponds to the situation where two pho- 
tons with different frequencies loq and u>\ are incident 
into the Raman medium. The input state of the fields 
can be written as 

l*in> = |l)o|l)a II |0>,=&S(0)St(0)|0>. (45) 

8^0,1 




FIG. 7: Joint probability Wo\ for finding one photon in each 
of the sidebands and 1 as a function of the effective medium 
length k,L. 

The output state of the fields is given by |\&out) = 
bl(-L/c)b\(-L/c)\0). With the help of Eq. Q, we find 

|tf out ) = -i^Y^e 1 ^ J q {nL)J q ^{nL)\2 q ) 

Q 

-iJ2^ (k+l)w/2 [MKL)Ji-i( K L) 

k<l 

+ J l (KL)J k _ 1 ( K L)}\l k l l ). (46) 

Here, the Fock state \2 q ) is the state of two photons in 
the sideband q with no photons in the other sidebands, 
and the Fock state |lfcl;) is the state in which there is one 
photon in each of the sidebands k and I but no photons 
in the other sidebands. 

It follows from Eq. I|46|) that the probability for finding 
two photons in the sideband q is 

= 2J 2 ( K L)J 2 _ 1 (kL). (47) 

The joint probability for finding one photon in each of 
the sidebands k and I (k ^ I) is given by 

W kl = [J fc (KL)J,_i(/eL) + Ji{kL) J fc _ 1 («L)] 2 . (48) 

The probability W q = J2i^ q W q i for having one and 
only one photon in the sideband q is 

= J 2 (kL) + J^(kL) - 4J 9 2 ( K L)J 8 2 _ x ( Kj L). (49) 

The mean photon number of the sideband q is 

{h q )=J 2 q {KL) + J 2 q _ x {KL). (50) 

We find the relations W (2) q = W$ q , = W$ q , and 

(fi-q) = (hi +q ), which reflect the symmetry of the gen- 
erated Stokes and anti-Stokes sidebands with respect to 
the two input sidebands and 1. 

When we insert k = and I — 1 into Eq. ^48(1. we 
obtain the following expression for the joint probability 
for finding one photon in each of the sidebands and 1: 

Woi = [J$(kL)-J 2 ( K L)] 2 . (51) 
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This expression shows that the joint probability Wqi may 
become zero at certain values of kL [see Fig. [7j . This is 
a signature of destructive interference between two chan- 
nels that form the state | lo 1 1 ) ■ In the first channel, each 
of the two photons individually transmits through the 
medium without any changes. The two-photon prob- 
ability amplitude for this channel is Jq(kL)Jo(kL) = 
Jq(kL). In the second channel, both the photons are 
scattered from the prepared Raman coherence and ex- 
change their sidebands. The two-photon probability 
amplitude for this channel is e l7T ^ 2 Ji(KL)e" r / 2 J\{kL) — 
— Ji(kL). Since Raman scattering produces a phase shift 
of 7r/2 for each photon, the probability amplitudes for the 
two channels (the transmission and scattering of both 
the photons) are 180° out of phase. The interference be- 
tween the two channels is therefore destructive, yielding 
the output state |loli) with the joint probability Woi 
given above. When the medium length L is such that 
Jq(kL) — Ji(kL), the interference between the two two- 
photon amplitudes becomes completely destructive, and 
therefore the state | IqIi) is removed from the output 
state (@I)J). We denote such a medium length by Lf. The 
positions of the zeros of Wqi depicted in Fig. [7| indicate 
that the first three values of Lf are given by nLt — 1.44, 
3.11, and 4.68. It is interesting to note that kL/ can 
be determined in an experiment using a single-mode in- 
put. Indeed, in the case where a single sideband is ini- 
tially excited, the mean photon numbers of the generated 
sidebands are given by Eq. (JSJ. Therefore, the effective 
medium length kL t corresponds to the situation where 
the probe sideband and its adjacent sidebands ±1 have 
the same mean photon numbers at the output. 

There exist literatures on two-photon interference in 
various systems |0, 0, |23|| . Two-photon interference in 
coherent Raman scattering, described above, is an anal- 
ogy of two-photon interference at a conventional beam 
splitter 0, [l^. We emphasize that two-photon inter- 
ference in coherent Raman scattering involves copropa- 
gating photons with different frequencies in a collinear 
scheme. 



D. General relation between the P representations 
of the input and output states 

To be more general, we consider the case where an ar- 
bitrary number of sidebands is initially excited. We find 
that an arbitrary multimode coherent state |{«g(0)}) of 
the input fields produces a coherent state \{a q (L/c)}) 
of the output fields. Here, the output amplitudes 
{a q (L/c)} are linearly transformed from the input am- 
plitudes {a g (0)} as given by 

a q (L/c) =J2e t{q - q ' )n/2 J q ^(nL)a q/ (0). (52) 
i' 

Consequently, the diagonal coherent-state representation 
P\rA\ctqY) °f an arbitrary input quantum state pi n deter- 
mines the representation P nt({o: q }) of the output state 



Pout via the equation 

Pout ({« J) = Pn(K}). (53) 
Here, we have introduced the notation 

If the input state p- m is a classical state [l^, Pi n ({a q }) 
must be non-negative and less singular than a S func- 
tion, and consequently so must P ou t({ag})- In this case, 
the output state /5 ou t is also a classical state. Moreover, 
since the multimode coherent state is separable 

and the weight factor P ut({o-q}) is non-n egative, the out- 
put state y0 ou t is, by definition, separable |24j,|25j. There- 
fore, a necessary condition for the output fields to be in 
an inseparable (entangled) state or, more generally, in a 
nonclassical state is that the input field state is a non- 
classical state. A similar condition has been derived for 
the beam splitter entangler [l7|. Note that, in the case 
where we use a single-mode input field q — 0, prepared in 
an arbitrary quantum state with the coherent-state rep- 
resentation Pj^ (a) (the Stokes and anti-Stokes sideband 
fields are initially in the vacuum state), Eq. becomes 

P out ({a q }) = P^(a' )l[S(a' q ). (55) 

It has been shown in Ref. that, when the input field is 
prepared in an even or odd coherent state, a multipartite 
entangled coherent state can be generated. 



V. CONCLUSIONS AND DISCUSSIONS 

We have studied the quantum properties of multi- 
order sidebands generated by the beating of a quantum 
probe field with a prepared Raman coherence in a far- 
off-resonance medium. Under the conditions of negligi- 
ble dispersion and limited bandwidth, we have derived a 
Bessel-function solution for the sideband field operators. 
We have analytically and numerically calculated various 
quantum statistical characteristics of the multiorder side- 
band fields. 

We have examined the quantum properties of the side- 
band fields in the case of single- mode quantum input. Wc 
have shown that, when we change the effective medium 
length or the Raman sideband order, the autocorrelation 
functions, the cross-correlation functions, the photon- 
number distributions, and the squeezing factors undergo 
oscillations governed by the Bessel functions. When the 
sideband order is higher, the onset of the sideband gen- 
eration occurs later and, therefore, so does the onset of 
the sideband autocorrelation functions. The mean pho- 
ton number and the autocorrelation functions of each 
sideband reach their largest values at the same optimal 
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medium length determined by the first peak of the cor- 
responding Bessel function. The higher the sideband or- 
der, the larger is the optimal length and the smaller is 
the maximal output values of the mean photon number 
and sideband autocorrelation functions. 

Meanwhile, the normalized autocorrelation functions 
and normalized squeezing factors of the probe field are 
not altered by the parametric beating process. They are 
replicated into the comb of generated multiorder side- 
bands. As the result, the single-mode input field and 
the generated sidebands have identical normalized auto- 
correlation functions and identical normalized squeezing 
factors. In other words, they have similar quantum sta- 
tistical properties - the same type of photon statistics 
and the same type of squeezing. In addition to this re- 
semblance, it has been shown that an input field in a 
coherent state can produce sideband fields in coherent 
states. It has also been shown that, when the input field 
is prepared in a thermal state, the reduced state of each 
generated sideband is also a thermal state. Therefore, 
the multiorder coherent Raman process can be used to 
multiplex the statistical properties of a quantum probe 
field into a broad comb of different frequencies. 

As far as replicating the statistical properties of the 
input probe into its sidebands is concerned, the Raman 
medium appears to behave as a linear system. However, 
the replication of the normalized autocorrelation func- 
tions and normalized squeezing factors of the probe field 
does not mean the replication of the quantum state at 
all. The photon-number distributions and the quantum 
states of the sidebands evolve in a rather complicated 
way. Cross-correlations between the sidebands can be 
generated from initially uncorrelated fields. An insepa- 
rable state can be generated from a separable nonclassical 
state. Although the dynamics of our model system is lin- 
ear with respect to the field variables, the possibilities of 
interesting quantum phenomena such as anti-correlation 
generation, squeezing multiplexing, and entangled-state 
generation represent the quantum properties that cannot 
be described by the classical statistics of the fields with 
positive P functions. 

We have also studied the mixing and modulation of 
photon statistical properties in the case of two-mode 
quantum input. We have shown that the prepared Ra- 
man coherence and the medium length can be used as 
control parameters to switch a sideband field from one 
type of photon statistics to another type, or from a non- 
squeezed state to a squeezed state and vice versa. In 
addition, we can switch nonclassical properties, such as 



sub-Poissonian photon statistics and squeezing, from one 
frequency to another frequency. We have also shown an 
example of quantum interference between the probability 
amplitudes for a pair of photons with different frequen- 
cies. 

We have made interesting observations that the mul- 
tiorder coherent Raman scattering behaves in many as- 
pects as a conventional beam splitter and hence can be 
called a multiorder Raman beam splitter. The two sys- 
tems have the same underlying physics: the fields are 
linearly transformed from the input values. However, 
the two systems are different in their natures. Unlike 
the conventional beam splitter, the multiorder coherent 
Raman process can efficiently produce a broad comb of 
sideband fields whose frequencies are different and are 
separated by integer multiples of the Raman modulation 
frequency. The number of generated Raman sidebands 
increases with the effective medium length, which is pro- 
portional to the product of the medium length and the 
prepared Raman coherence. The Bessel functions of the 
effective medium length play a similar role as the trans- 
mission and reflection coefficients of a conventional beam 
splitter. 

The ability of the far-off-resonance Raman medium to 
generate a broad comb of fields with similar quantum 
statistical properties and to switch the quantum statis- 
tical characteristics of the radiation fields from one type 
to another type may find useful applications for high- 
performance optical communication networks. In addi- 
tion, two-photon interference in coherent Raman scat- 
tering may find various applications for high-precision 
measurements and also for quantum computation. 

Finally, we emphasize that the coupling between the 
Raman sidebands can be controlled by the magnitude of 
the prepared Raman coherence, that is, by the intensi- 
ties of the driving fields. In a realistic far-off-resonance 
Raman medium, such as molecular hydrogen or deu- 
terium vapor 0, |p , solid hydrogen 0, |g , and rare-earth 
doped dielectrics [Hj, a large Raman coherence and, con- 
sequently, a large number of Raman sidebands can be 
generated by the two-color adiabatic pumping technique. 
In such a system, the generation of a broad comb of 
high-order Raman sidebands with nonclassical statisti- 
cal properties is, in principle, feasible. Therefore, we 
expect that the coherent Raman scattering technique us- 
ing quantum fields will become a practical and efficient 
method for a wide range of applications in nonlinear and 
quantum optics. 
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